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JT) (57) Abstract: The invention concerns a method for improving depth discrimination in optical reproduction systems, said method 
being particularly useful in optical microscopy for enhancing the image quality when examining three-dimensional elongated ob- 
qq jects. The invention concerns in particular the method of structured illumination such as disclosed in WO 97/6509. The invention 
ff^ is characterized in that it consists in determining the influence exerted by the variations of luminosity of the light source and the 
discoloration of the object in the case of a fluorescent illumination, said factors being taken into consideration when calculating the 
structure of the object. 



(57) Zusammenfassung: Die Erfindung betrifft ein Verfahren zurTiefendiskriminierung optisch abbildender Systeme. Es ist insbe- 
sondere in der Lichtmikroskopie anwendbar, um eine Erhohung der Bildqualitat bei der Untersuchung dreidimensional ausgedehnter 
Q Objekte zu erreichen. Sie ist anwendbar beim Verfahren der strukturierten Beleuchtung wie es in WO 97/6509 beschrieben ist. Dazu 
£^ werden Einfliisse durch Helligkeitsschwankungen der Lichtquelle, der Positionierung der abgebildeten periodischen Struktur und 
^ Ausbleichen des Objekts bei Fluoreszenzbeleuchtung ermittelt und bei der Berechnung der Objektstruktur beriicksichtigt. 
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Verfahren zur Verbesserung der Tiefendiskriminierung optisch abbildender Systeme 

Die Erfindung betrifft ein Verfahren zur Tiefendiskriminierung optisch abbildender Systeme. 
Es ist insbesondere in der Lichtmikroskopie anwendbar um eine Erhohung der Bildqualitat 
5 bei der Untersuchung dreidimensional ausgedehnter Objekte zu erreichen. 

In der klassischen Lichtmikroskopie besteht bei der Untersuchung dreidimensional 
ausgedehnter Objekte, d.h. von Objekten, deren Ausdehnung entlang der optischen Achse 
grower ist als die Tiefenscharfe der vsrwendeten Objektive, das Problem, dass das scharfe 
10 Bild aus dem Fokusbereich mit auBerfokalen Bildanteilen, welche demzufolge unscharf 
abgebildet werden, uberlagert ist. 

Um dieses Problem zu iiberwinden ist die Nutzung der konfokalen Abbildung bekannt, bei 
der mittels eines sogenannten Pinholes das aus dem nicht im Fokus liegenden Bereich 
stammende Licht ausgeblendet wird und somit nicht zur Abbildung beitragt Auf diese Weise 
1 5 entsteht ein sogenannter optischer Schnitt. Diese konfokale punktweise Abbildung erfordert 
das Abrastern des Objektes in der Bildebene um ein Bild zu erhalten. Diese Abtastung kann 
entweder mittels Scannern in Laser Scanning Mikroskopen oder mittels Nipkow-Scheiben 
erfolgen. 

Durch Aufiiahme mehrerer optischer Schnittbilder in verschiedenen Fokuslagen kann ein „z- 
20 Stapel" gewonnen werden, welcher eine dreidimensionale Darstellung des Objektes 
ennoglicht 

Ein anderer Weg zur Erzeugung optischer Schnitte ist die Anwendung der strukturierten 
Beleuchtung. Diese wurde erstmals yon Meir Ben-Levy und Eyal Pelec in WO 97/6509 

25 angegeben. Verbesserungen und Erweiterungen dieser Methode wurden von Tony Wilson et 
al. in WO 98/45745 und von Volker Gerstner et al. in WO 02/12945 beschrieben. Auf die 
Offenbarung dieser drei Schriften wird hiermit als Referenz verwiesen. 
hi WO 97/6509 wird das Objekt mit einer periodischen Struktur (Sinus- oder Rechteckgitter) 
beleuchtet und mittels einer Kamera ein Bild des Objektes aufgenommen, digitalisiert und in 

30 einem Speicher abgespeichert. AnschlieBend wird die periodische Struktur innerhalb der 

Bildebene so verschoben, dass die Phasenlage der Struktur geandert wird, und wieder ein Bild 
aufgenommen und abgespeichert. Dieser Vorgang (Verschiebung, Bildaufiaahme, 
Abspeicherung) kann mehifach wiederholt werden. Anschhei3end werden die Bilder 
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miteinander verrechnet urn ein Schnitfbild zu erzeugen. Als mathematischer Aasatz wird 
dabei eine Fourierentwicklung angegeben, was zu einem komplizierten Formelapparat fuhrt 
In WO 98/45745 wird eine einfachere Formel fiir die Schnittbilder angegeben, welche sich 
durch Vereinfachung des Ansatzes aus WO 97/6509 fur den Falls jeweils gleicher 
5 Phasenverschiebungen zwischen den einzelnen Aufhahmen herleiten lasst 

Die Realisierung der fur die Anwendung der angegebenen Verfahren notwendigen 
physikalischen Randbedingungen erweist sich in der Praxis als sehr schwierig. So fuhrt z. B. 
die Schwankung der Lampenhelligkeit zwischen den verschiedenen Aufhahmen zu 
streifenformigen Artefakien in den generierten Schnittbildera. Bei Flucreszenzobj ekten treten 
1 0 zusatzliche Probleme durch die zeitabhangige Ausbleichung der Fluoreszenzfarbstoffe auf, 
welche gleichfalls zu Fehlern fuhrt. Auch die notwendige Konstanz der einzelnen 
Phasenverschiebungsschritte lasst sich in der Praxis nicht einhalten. 
In WO 02/12945 wird daher vorgeschlagen, mittels Auskopplung eines Teils des zur 
Beleuchtung des Objekts dienenden Lichtes, Registrierung der Intensitat und anschlieflende 
1 5 Normierung der Einzelaufiiahmen den Einfluss der zeitlich unterschiedlichen 
Lampenhelligkeit auszugleichen. Fiir die Beriicksichtigung ungleicher 
Phasenverschiebungsschritte wird ein Gleichungssystem (dort Gl. 22) angegeben. Zum 
Ausgleich der Ausbleichung wird vorgeschlagen, statt der mindestens notwendigen drei 
Aufhahmenje Schnittbild 6 Aufhahmen in der Reihenfolge 1-2-3-3-2-1 zu 
* 20 registrieren, zwei Schnittbilder (aus 1-2-3 und 3 - 2 - 1) zu berechnen und von diesen 
beiden den Mittelwert zu bestimmen. 

Die Realisierung dieser Vorschlage erfordert erheblichen instrumentellen Aufwand. 
AuiJerdem verlangert die Aufiiahme zusatzlicher Bilder die notwendige Aufhahmezeit und 
erhoht damit auch die Belastung der Probe durch Beleuchtung mit dem 

25 Fluoreszenzanregungslicht 

Die Erfindung stellt sich die Aufgabe, die Nachteile des Standes der Technik zu vermeiden 
und ein verbessertes Verfahren zur Bestknmung optischer Schnittbilder anzugeben. 
. Diese Aufgabe wird durch ein Verfahren zur Verbesserung der TiefendiskrimMerung optisch 
abbildender Systeme gemaB dem Hauptanspruch gelost. Vorteilhafle Ausgestaltungen sind in 

30 den Unteranspriichen angegeben. 

Die Erfindung wird im folgenden anhand der Figuren naher beschrieben. 
Es zeigen: 

Fig. 1 den prinzipiellen Aufbau eines Mikroskops mit strukturierter Beleuchtung 
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Fig. 2 eiae schematische Darstellung der erfindungsgemaflen Bandpassfilter in der 
Fourierebene 

Fig. 3 schematisch die Lage der Abbildungen der periodischen Struktur 

Fig. 4 einen Ablaufplan fur die Bestimmung der Koirekturwerte durch lineare Optimierung. 

5 

In Fig.l ist das vereinfachte optische Schema der sfrukturierten Beleuchtung am Beispiel der 
Durchlichtbeleuchtung dargestellt Eingezeichnet ist der bildgebende Strahlengang 
(Lukenstrahlengang). Eine 1-dimensionaleperiodische Struktur (Transmissionsgitter) (3), 
welche sich in eiaer Scharfeebene der dargestellten cptischen Anordnung befindet, wird durch 

10 eine Lichtquelle (1) und nachfolgender Kollektoroptik (2) beleuchtet. Dem Gitter folgt in 
Lichtrichtung eine planparallele Glasplatte (4). Der Winkel der planparallelen Platte zur 
optischen Achse kann definiert eingestellt werden. Durch die nachfolgende 
beleuchtungsseitigen Optiken (5, 6) (Kondeaisor) wird die Struktur in die Praparatebene (7) . 
abgebildet. Von dem Praparat ausgehendes Licht wird wiederum uber ein Linsenpaar (8, 9) 

1 5 (Objektiv und Tubuslinse) in die darauffolgende Scharfenebene (10) abgebildet, in welcher z. 
B. eine CCD-Matrix einer digitalen Kamera angeordnet sein kann. Dabei dient die 
planparallele Glasplatte (4) durch definierte Verkippung zur Verschiebung der Abbildung der 
Gitterstruktur (3) auf das in der Praparatebene (7) befindliche Objekt. Vorzugsweise dient im 
Fall der Auflicht-Fluoreszenz-Beobachtung das Objektiv (8) gleicbzeitig als Kondensor. 

20 Wie in WO 02/12945 beschrieben wurde werden jetzt mit mindestens 3 Stellungen der 

Glasplatte (4) mittels der digitalen Kamera jeweils eine Registrierung der sich ergebenden 
Helligkeitsverteilung im Objekt durchgefuhrt. 

Dabei kann die Helligkeitsverteilung Ij (x, y) for den Fall eines Sinus/Kosinus-Gitters 
vereinfacht durch die Formel 

25 

7 l (x,y)=/ 0 (x 5 j;)<l+m(x,y).cos(^ 0 (x > j;)+^(x,j;)) (1) 

angegeben werden, wobei i=0 ... N-l die i-te Phasenposition des projizierten Gitters und N die 
Anzahl der Registrierungen ist, m(x, y) die Modulationstiefe des Objekts (und damit die 
30 gesuchte Bildinformation an der Stelle x, y) ist und die <fo die Phasehwerte sind. 

Diese Gleichung enthalt die drei unbekannten GroJJen Io, m und <J>q. Somit Vann durch 
mindestens drei Messungen bei gezielt variiertem §\ (i = 1, 2, 3) die Bestimmung dieser 
GroBen erfolgen. Die Losung kann durch einen kleinste-Quadrate-Ansatz aus den Messimgen 
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gewonnen werden. Dazu bietet es sich an, Gleichung (1) in einer kompakteren Form 
darzustellen und die Cosinus-Funktion mittels Additionstheorem umzuschreiben: 

/(x.y,A) = « 0 (*.y)+«i(*.J')'/i(A) + «»(*.J')-/a(*) ( 2 ) 

5 

mit 



/i (A ) = «»#, 
/i(^,) = sin^ 
a 0 (.x,y)=I Q (x,y) 

«i (x, y) = I 6 {x, y) • m(x, y) ■ cos <f> 0 (x, y) 
a 2 (x, y) = -I 0 (x, y) • m(x, y) ■ sin <f> 0 (x, y) 

10 Die Funktionen fj und hangen somit nur von den im Prinzip frei zu wahlenden 

Phasenverschiebungen i?\ ab. Fur den Fall nichtsinusformiger, aber periodischer Struktur kann 
die Helligkeitsverteilung Ii (x, y) auch durch Reihenentwicklung approximiert werden, das 
Prinzip der Rechnung bleibt dabei im Grunde gleich. Die kleinste-Quadrate-Losung lautet in 
Matrixschreibweise 

15 

M>a=b (3) 
nut 



20 M = 



(4) 



V N 



mit N = Anzahl Messungen (bier: Phasenschritte) und 



«i(*>J') 



(5) 
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sowie 



b = 



(6) 



5 Die Losung dieser Matrixgleichung wird mathematisch durch invertieren der Matrix M 
gewonnen. 



(7) 



10 



Somit konnen aus den Messungen, dargestellt durch den Vektor b, und den 
Phasenverschiebungen des Gitters die gesucbten unbekannten Grofien ennittelt werden. 
Ihsbesondere ergibt sich die Modulationstiefe am Ort (x, y) zu 



15 



20 



m(x,y) = 



<*o( x >y) 



(8) 



Da mit Defokussierung des Objektes der Gitterkontrast und damit zwangslaufig die 
Modulationstiefe im Objekt abnimmt, stellt Gleichung (8) eine Bestimmungsgleichung zum 
Erhalt eines tiefendiskriminierten Bildes (optischer Schnitt) dar. Gleichungen (2) und (3) 
liefern einen universalen Losungsweg, der die Anpassung der Phasenlagen fa an die jeweils 
experimentellen Gegebenheiten erlaubt So ist etwa die Wahl der Zahl der Schritte und die 
Positionen im Objektraum (Ort der abgebildeten Gitterlinien) frei wablbar. 



Jedes optische System ubertragt gernafi seiner optischen Transferfimktion Raumfrequenzen 
bis zu einem Grenzwert, der duroh die numerische Apertur NA der Optik und der 
25 Lichtwellenlange X bestimmt ist: 



k = 



4nNA 
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Dies hat zur Folge, dass von einem in den Objektraum abgebildeten Gitter, welches kein 
reines Sinus-Gitter ist, die hoheren Hannonischen der Gitterperiode bis zu eben dieser 
Grenzfrequenz lw iibeitragen werden und sich in dem Fourier-transfonnierten Bild des 
5 Objektes als lokale Maxima wiederfinden. Auch in einem nach dem oben angefuhrten 

Schema ermittelten Schnittbild sind diese Komponenten des Gitters enthalten und resultieren 
in Streifenartefakten. Diese Artefakte sind dadurch eliminierbar, dass die Hannonischen des 
projizierten Gitters durch Fourier-Transformation im Raumfrequenzspektrum lokalisiert imd 
durch gezieltes Setzen von Filtern (Bandpass-Filter, Bandstops) entfernt werden. Eine 

10 anschlieBende Riicktransformation fiihrt dann zu artefaktfreie Bildem. Diese Mefhode ist 
schematisch in Fig. 2 dargestellt hi der Fourierebene (1 1) sind mit (12) die hoheren 
Hannonischen der Gitterstruktur (3) bezeichnet, (13) bezeichnet den Nullpunkt und entspricht 
dem Gleichanteil der Beleuchtung, d.h. einer gleichmaBigen, nicht strukturierten Beleuchtung. 
Durch Einbringen den Hannonischen (12) entsprechender Bandpass-Filter in den 

1 5 Strahlengang tragen diese Bereiche nicht zur Abbildung bei und werden damit nicht von der 
digitalen Kamera registriert. 

Bei fluoreszierenden Objekten tritt ein weiterer Grund zum Auftreten von Streifenartefakten 
hervor, namlich das Fluoreszenz-Bleichen (photobleaching, fluorescence fading). Hierbei 
kommt es durch die Abbildung der Struktur (z.B. des Gitters) in das Objekt zu Regionen 

20 unterschiedlich starker Ausbleichung der Fluoreszenz und somit schlieBlich zu 

Streifenartefakten in jedem in der Folge aufgenommenem Bild. Die Frequenzanteile der 
Artefakte im Spektrum eines solchen Bildes sind besonders stark bei der Grundfrequenz der 
projizierten Struktur. Durch Entfernen der Raum frequenzanteile der projizierten Struktur im 
Fourier-transformierten Bild durch entsprechende Bandpass-Filter und anschlieBender 

25 Rucktransformation kann auch hier prinzipiell ein artefaktfreies Bild erzeugt werden. 

Die Verhinderung von diesen bei der Fluoreszenzbeleuchtung auftretenden Artefakten lcann 
auch durch die Bestimmung und Ausnutzung der Bleicheigenschaften von fluoreszierenden 
Objekten bei der sequentiellen Aufiiahme von Bildern desselben Objekts erreicht werden. 
30 Hierbei ergibt sich pro Bildaufhahme durch die Beleuchtung eine Bleichung des Objekts an 
den beleuchteten Stellen. Die Starke des Bleichens hangt neben der Beleuchtungsintensitat 
von den photochemischen xind photophysikahschen Eigenschaften des fluoreszierenden 
Objektes ab. Die Bestimmung der Fluoreszenzstarke des Objektes durch Integration der 
Intensitaten uber alle oder ausgewahlte Objektorte (x, y) ennoglicht die ReskaUerung der 
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Bildsequenz nach erfolgter Au f hahme der Bilder. Dazu werden entweder iiber analoger 
Kamera und Analog-Digital- Wandler oder einer CCD-Kamera in einem Computer 
gespeicherten Bildinfonnationen durch paarweise Ermittlung von Quotienten der 
Fluoreszenzstarken der Bilder die Intensitaten an den Bildorten (x,y) auf einen Maximalwert 
5 normiert. Somit lasst sich der Fluoreszenzverlust durch Fluoreszenzbleichen berucksichtigen. 
Dazu konnen nach Fig. 3 Regions of Interest (ROI) definiert werden, die bei beiden 
Bildaufhahmen gleichennafien beleuchtet werden und in Folge fluoreszieren, das MaB der 
Abnahme der Fluoreszenz in diesen Bereichen ist ein MaB fur die Bleich-Rate und kann damit 
ffir die Normierung der verschiedenen Registrierungen der Heliigkeitsverteilung 

10 herangezogen werden. Die Definition dieser ROIs kann dabei sowohl vom Benutzer iiber 

entsprechende Eingabemittel als auch automatisch (aus der jeweiligen Lage der abgebildeten 
Struktur) erfolgen. In Fig. 3 stellt (14) die Verhaltnisse bei einer ersten Registrierung der 
Heliigkeitsverteilung auf dem Objekt dar, (15) die Verhaltnisse bei einer zweiten 
Registrierung mit veranderter Phasenlage der abgebildeten Struktur (16). Die Regions of 

1 5 Interest (17) werden vorteilhaflerweise so definiert, dass in beiden Registrierungen 
gleichmafiige Beleuchtung erfolgt. 

Das erfindungsgemafie Verfahren erlaubt die Bestimmung von optischen Schnittbildern 
praktisch fur beliebig gewahlte Phasenverschiebungen fa. Zur Verhinderung von Fehlern sind 
20 die Phasenverschiebungen <j>i mit hoher Genauigkeit zu bestimmen. Dieses lasst sich in 

einfacher Weise mit dem im folgenden beschriebenen Kalibrierungsverfahren sichern, bei 
dem die Einflusse des gesamten Systems beriicksichtigt werden. 

Dazu wird zunachst ein Referenzbild (Objekt mit iiberlagertem Gitter) aufgenommen und 
registriert. Daraufhin erfolgt eine geringe Phasenverschiebung des Gitters auf eine zunachst 

25 unbekannte Phasenlage des Gitters (es soli diese „Verschiebung pro Ansteuerungssignal" erst 
kalibriert werden) und erneute BildaujBoahme. Die erhaltenen zwei Bilder werden nun 
rechnerisch verglichen. Dies kann ohne Einschrankung der Allgemeinheit durch 
Differenzbildung, Summierung oder jede andere Art erfolgen, die auf eine Merit-Funktion 
fuhrt. Resultat ist zunachst im Allgemeinen ein Streifen-behaftetes Bild. Die Schritte: 

30 . Verschiebung des Gitters um einen kleinen Betrag, Bildaufeahme, Vergleich mit dem 

Referenzbild werden so lange wiederholt, bis die Merit-Funktion den Extremwert erreicht, 
also z.B. ein Differenzbild nur noch aus.Untergrundrauschen besteht oder ein Summenbild 
maximale Werte erreicht Ist letzteres erfullt, ist die KaUbration abgeschlossen und ein so 
erhaltener Ansteuerwert kann in einem Speichermedium zum spateren Zweck abrufbar 
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gehalten werden. Dieser Wert entspricht damit genau der Verschiebung des Gitters um eine 
voile Periode. Alternativ ist es auch moglich, ein Summenbild auf gleichmafiige 
Helligkeitsverteilung (dJi. Verschwinden der Streifenstruktur) zu bewerten. In diesem Fall 
entspricht dieser Wert dam der Verschiebung des Gitters tun eine halbe Periode. 
5 Vorzugsweise wird diese Prozedur mit einem Spiegel als Objekt durchgefiihrt Der 
beschriebene Vorgang kann sowohl manuell als auch automatisiert erfolgen. 

. Im folgenden wird an Hand der Fig. 4 ein alternatives Verfahren zur Bestimmung der 
Einfiiisse von ktaipenhemgkeitsschw Ausbleichen des Objekies bei 

10 Fluoreszenzbeleuchtung und nicht-sinusfonnigem Gitter beschrieben. 

Dazu wird die zeitliche Veranderung der Beleuchtung durch intensitatsproportionale skalare 
Faktoren 9 t > 0 beschrieben werden, welche bei jeder aufeinanderfolgenden Aufioahme 
beriicksichtigt werden (die Abhangigkeit der Grofien g i? o u usw. vom Ort (x,y) wird im 
1 5 Folgenden zur Vereinfachung nicht mehr explizit angegeben) : 

gi=0A (9) 
Die modellierte Beobachtung & ergibt sich nun aus dem Produkt mit dem idealen, 
konsistenten Phasenbild o im Um vollstandig ein gestortes System zu beschreiben, braucht man 
rxur N-l Faktoren Q t zu bestimmen. 

20 

Das Ausbleichen des Objekts bei Fluoreszenzbeleuchtung wirkt in Abhangigkeit von der 
einwirkenden Strahlungsintensitat uber die Zeit Diese ist wegen des projizierten Gitters 
ortlich verschieden und erfordert deshalb eine Behandlxmg, die die Gitterfiinktion 
beriicksichtigt. Eine einfache Beschreibung wird moglich, wenn man das Ausbleichen als eine 
25 lineare Funktion der Beleuchtungsintensitat annimmt: 

k^fl-rfo^f^^k^^l mitk 0 =l (10) 

Der aus ortsvarianten Faktoren bestehende Vektor 0 < k/ < 1 soil abhangig vom Verlauf der 
Bestrahlungsintensitat und einem wahlbaren Mafi 1 > d > 0 fur das Ausbleichen die 
Abschwacbung beschreiben wenn man ihn komponentenweise mit den jeweils noch nicht 
30 ausgeblichenen Beobachtungen 0/ multipliziert: 

8/=0/k / (11) 
Aus (10) geht hervor, dass immer jede Beobachtung mit den Faktoren aller vorhergehenden 
Phasenaufiiahmen multipliziert wird, diese sind bereits z.T. ausgebhchen und uberlappen sich 
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mit der momentanen Phase. Altemativ kann man k, auch exponentieU d.h entsprecheaid der 
. Abklingkurve eines bestimmten Fluorochromes beschreiben. 

• Nicbtsinusformige Gitterverlaufe verursachen Streifenartefakte insbesondere fur ganzaaMige 
5 Vielfeche der Gitterfrequenz. Unter genauer Kenntnis von Frequenz, Phasenlage und 

Kurvenform des Gitters bei den einzelnen Aufhahmen wird es moglich einen sinusformigen 
Verlauf zu erzeugen. Ein Korrekturvektor 1^ ahnlich wie unter (10), hier jedoch entlang der 
Gitteiperiodizitat, soil unter (12) betrachtet werden. 

f. 

K °" < 12 > 
10 Die Nachbildung der observierten Gitterfunktion ist mit/^ gekennzeichnet. Ein ideal 

sinusformiger mtensitatsverlauf wird mit/^ zum Ausdruck gebracht. Die Konstante s im 
Nenner soli die Division durch sehr kleine Werte von f Gltter ubemehmen. Die Nachbildung der 
observierten Gitterfunktion kann auf verscbiedenen Wegen geschehen. Eine einfache 
Moglicbkeit einer Approximation bietet sich z.B. durch die Verwendung einer Fourierreihe 
15 fur die Trapezfunktion an. 

Mit b kann man die Flankensteilheit der Trapeziibergange variieren, x ist die Ortskoordinate 
und M die Anzahl der benutzten Fourierkoeffizienten. Altemativ lassen sich auch genauere 
ModeUe verwenden. Diese erfordem die Faltung der wahren Gitterfunktion mit einer lateralen 
20 Komponente der PSF des Mikroskops. Von entscheidender Bedeutung ist eine phasen- und ' 
frequenzsynchrone Synthetisierung des Korrekturvektors \ Kon . Ist dies gewahrleistet, so kann 
ahnlich wie unter (11) mit einer einfachen komponentenweisen Multiplikation ein korrigiertes 
Ergebnis fur die Observation erzielt werden. 

%Korr = 8 l^iXorr (14) 

25 DiekorrigiertenPhasenbildererhaltmannunaus: 

h =°, = g !*'% rr , mit 9, > 0 n k,-> 0 . (15) 

Durch Einsetzen von (15) in Gleichung (3) bzw. (8) erhalt man eine allgememe Formel zur 
Rekonstruktion, welche auBerdem die Parameter 9 t , d und b enthalt 
30 Die Losung des gesteUten Problems kann beispielsweise durch Umformung in eine 

Extremalaufgabe erfolgen. Mit den Hilfsmitteln numerischer Optimierung kann man die 
Parameter 0, , d und b unter Zuhilfenahme einer Merit-Funktion bestimmen. 
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M(0 f ;rf;6) = a o |F{5} o |>^ 

Hierbei ist F { a } co eine Komponente der Fourier-Transformation des Resultatsvektors a . 
. aus Gleichung (5), Dabei sind neben der Fourier-Transformation auch andere 
5 Funktionaltransfonnationen geeignet. 

Dabei ist die Meritfunktion (Gleichuag (16)) hier nur beispielsweise genannt, es sind auch 
andere Meritfunktionen anwendbar. 

Ziel ist es 0, 9 d und b so- zu variiereu, dass (16) minimiert wird. Dafiir bietet sich eine 
Vielzahl von Methoden an die keinesfalls eingeschrankt zu Gradienten- oder 
10 Liniensuchmethoden zu sehen sind. Jeder TransfonnationskoefiBzient ist auBerdem gewichtet 
mit einem Koeffizienten «/. Darait kann der Algorithmus an verschiedene Signal-Rausch 
Abstande oder bevorzugte Frequenzen angepasst werden. Bin diesbezuglich vorteilhafter 
Wert den Gleichanteil zu wichten ergibt sich mit folgender Beziehung: 

15 |ywj^g>j*-^g>.r m 

Die allgemeine Funktionsweise des hier beschriebenen Korrekturapparates ist dem 
schematischen Ablaufplan in Fig. 4 zu entnehmen. Die Verwendung von Projektionen ist 
vorteilhaft um effizient die Parameter 0, , d und b zu bestimmen. Die Parameterennittlung ist 
aber nicht auf diesen eindimensionalen Fall eingeschrankt, sondein funktioniert in der Ebene 

20 in analoger Weise. 

Es werden bei fester Fokuslage von der Kamera mindestens drei Helligkeitsverteilungen (1 8), 
(19), (20), ggf. auch mehr (21), mit verschiedenen Phasenlagen i=l, 2, 3, . . N der auf das 
Objekt abgebildeten Gitterslruktur registriert. Mittels Projektion werden diese 2- 
dimensionalen HeUigkeitsverteilungen jeweils in eine 1-dimensionale Verteilung (22), (23), 

25 (24), und ggf. (25) uberfuhrt. Mit Anfangswerten fur die Parameter 0, , d und b werden nach 
Formel (15) korrigierte Phasenbilder berechnet und aus diesen mittels des Ansatzes aus 
Foimel (3) ein erstes Bild des optischen Schnitts errechnet (26). Die Anfangswerte der 
Parameter 0, , tfund b konnen dabei vom Bediener eingegeben oder auch automatisch 
geschatzt werden. Fur diese Schatzung konnen bekannten Tools der Minimierung / 

30 Maximierung einer Merit-Funktion benutzt werden (z. B. Simplex, Marquardt, o.a.). Aus 

diesem Ergebnis wird nun in (27) die Merit-Funktion nach Formel (16) bestimmt und in (28) 
auf das Erreichen eines Minimums getestet. Ist das Minimum noch nicht erreicht wird eine 
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neue Schatzung. der Parameter 9 t , d und 6 durchgefiihrt (29) und mit (26) fortgesetzt. Fur den 
Fall des Erreichens des Minimums der Merit-Funktion, d.h. weitere Variation der Parameter 
9 i , d und b lassen ihren Wert wieder ansteigen, sind die optimalen Parameter bestimmt und 
in (30) wird das optische Schnittbild mittels der Formeln (15), (3) und (8) bestimmt und in 
5 (3 1) ausgegeben. AnschlieBend kann die Fokuslage verandert werden und nach dem 
dargestellten Verfahren ein weiteres optisches Schnittbild bestimmt werden usw.. Durch 
dieses Abrastem des Objektes in z-Richtung lasst sich ein z-Stapel von Schnittbildem und 
damit eine dreidimensionale Darstellung des Objekts generieren. 

Ein anderer Ansatz zur Berucksichtigung von Ausbleicherscheinungen des Objekts bei der 
1 0 Fluoreszenzbeleuchtung beruht auf der Bestimmung einer lokalen Kon^ekturfunktion. Dazu 
wird eine Korrekturfonktion Kt(x) 



[x+t/2 
_ Jx-f / 



= Zlln —— (18) 



angenommen, deren lokaler Wert uber eine Gitterperiode t = 2w / co gebildet wird. Der Ansatz 
dieser Fimktion kann auch die Abhangigkeit von y, der hier zur Gitteiperiode senkrecht 
1 5 angenommen Koordinate beriicksichtigen, was hier jedoch auch Grunden der 
Ubersichllichkeit weggelassen wird. Unter.Verwendung einer ortsvariablen 
Ausbleichfunktion Ofa) und durch geeignete Umfomungen mit der Naherung (siehe auch 
Fonnel (1) und (9)) g i (x) = 9 t (x)I 0 (x)[l + m cos(<a x + , laBt sich zeigen 

«/(*)= fX +rn l£k dassbeigleichmassiger 

0, (£)/ 0 (£) d% + m\^ n e t (£) J 0 (<?) cos(^ +n)dt 

20 Ausbleichfiinktion diese Koirekturfunktion proportional zu l/9j(x) wird. 

Es hat sich gezeigt, dass es von Vorteil ist, wenn die Mittelwertbildung iiber mehr als eine 
Gitterperiode ausgefuhrtwird. 

Weiterhin hat sich gezeigt, dass diese Korrekturfunktion in der Nahe von Kanten in der 
Ausbleichfonktion zu Uberschwingen neigt und damit zusatzliche Verfalschungen des 
25 Ergebnisses erzeugen kann. Um dem abzuhelfen hat sich bewahrt, das Auftreten des 

Uberscbwingens durch Vergleich mit einem Mittelwert zu erkennen und an diesen Stellen 
anstatt des verfalschten Wertes den aus den Formeln (9) und (16) vorher bestimmten globalen 
Wert l/9t einzusetzen. 
Die entsprechende Fonnel lautet dann: 
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K -^ = Mx) fur |*r,(x)-J^>c ( (£)^|<-* 



(19) 



l/0 t ansonsten 



Dabei ist es vorteilhaft, die Schwelle e als Prozentwert der Variation der lokalen 
Veranderlichkeit der Korrekturfunktion innerhalb des Bildes zu setzen, z.B. 5%. Anstatt von 
1/0/ konnen auch interpolierte Werte von k£x) 9 welche unterhalb der Schwelle im 
5 zweidimensionalen Raum berechnet wurden, verwendet werden. 

Dieses Verfahren beriicksichtigt somit ortsabhangig schwankende Ausbleicherscheinungen 
z.B. aufgrund verschiedener Fluorophore oder schwankender Eigenschaften selbiger. 
: Die Realisierung der Erfindung ist nicht an die dargestellten Ausfiihrungsbeispiele gebunden, 
durch fachroannische Weiterentwicklungen oder andere Ansatze zur Beriicksichtiguiig der 
10 Helligkeitsschwankungen der Lichtquelle oder des Ausbleichen des Objektes bei der 
Fluoreszenzbeleuchtung wird der Bereich der Erfindung nicht verlassen. 
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Aospruche 

1 . Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme bei dem 
a) ein Objekt mit einer periodischen Stxuktur beleuchtet wird, 

5 b) eine Registrierung der so entstehenden Helligkeitsverteilung erfolgt, 

c) die Phasenlage der periodischen Struktur verschoben wird, 

d) die Schritte a) bis c) sooft wiederholt werden bis mindestens drei 
Helligkeitsverteilungen registriert wurden, 

e) die fegisixierten Helligkeitsverteilungen initeinander verrechnet werden urn eine 
1 0 Objekfhelligkeitsvertettung zu erhalten, 

gekennzeichnet dadurch, dass 

f) die Verschiebungen der Phasenlagen aus Schritt c) registriert werden, 

g) Helligkeitsschwankungen der Beleuchtung erfasst werden 

h) Ausbleicherscheinungen des Objekts bei Fluoreszenzbeleuchtung bestimmt werden 
1 5 und die Verrechnung unter Beriicksichtigung der in den Schritten f) bis h) 

gewormenen Ergebnisse erfolgt. 

2. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach 
Anspruch 1, gekennzeichnet dadurch, dass Ausbleicherscheinungen des Objekts bei 

20 Fluoreszenzbeleuchtung dadurch bestimmt werden, dass aus vorzugsweise 

aufeinanderfolgend registrierten Helligkeitsverteilungen Gebiete gleicher 
Beleuchtungsintensitat bestimmt werden und die Helligkeit dieser Gebiete naiteinander 
verglichen wird 

25 3 . Verfahren zur Verbesserung der Tiefendisb±nimerung optischer Systeme nach 

Anspruch 2, gekennzeichnet dadurch, dass ein Normierungsfektor durch eine 
Quotientenbildung aus den registrierten Helligkeiten der Gebiete gleicher 
Beleuchtungsintensitat erfolgt und mit Hilfe dieses Nonnierungsfaktor die 
registrierten Helligkeitsverteilungen nonniert werden. 

30 

4. Verfahren zur Verbesserung der Tiefendiskrimiuierung optischer Systeme nach 
Anspruch 1, gekennzeichnet dadurch, dass die Bestimmung der 
Objekthelligkeitsverteilung m(x, y) durch die Gleichung 
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m{x,y) = 



10 



15 



erfolgt, wobei die Bestimmung des Vektors 



durch Losen eines Gleicbungssystems aus der Verknupfong der registrierten 
Helligkeitsverteilungen l{ (x, y, mit den registrierten Phasenverschiebungen <(>i 
erfolgt, wobei i fur die Anzahl der Registrierungen der Helligkeitsverteilung steht 

Verfahren zur Verbesserung der Tiefendiskrirninierung optischer Systeme nach 
Anspruch 3 und 4, gekennzeichnet dadurch, dass das Gleichungssystem die 
nonnierten Helligkeitsverteilungen mit den registrierten Phasenverschiebungen 
verkniipft 



6. Verfahren zur Verbesserung der Tiefendiskrirninierung optischer Systeme nach 

Anspruch 1, gekennzeichnet dadurch, dass eine Kalibrierung der Verschiebungen der 
Phasenlage erfolgt, wobei 

a) ein Objekt, welches vorzugsweise eine Spiegelflache ist, mit einer periodischen 
20 Struktur beleuchtet wird, 

b) eijie erste Helligkeitsverteilung registriert wird, 

c) die Phasenlage der periodischen Struktur urn einen bestimmten geringen Betrag 
verandert wird, welcher registriert wird , 

d) eine zweite Helligkeitsverteilung registriert wird, 

25 e) der Unterschied der beiden Helligkeitsverteilungen bestimmt und bzgl. 

streifenformiger Verteilung bewertet wird, 

f) die Schritte c) bis e) solange wiederholt werden, bis- die Bewertung des 
Unterschiedes zwischen erster und folgender Helligkeitsverteilung auf ein Extremum 
fuhrt,und 

30 g) der so gefundenen Wert der Phasenlage der periodischen Struktur registriert wird. 
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7. Verfahren zur Verbesserung der Tiefehdislaininiehing optischer Systeme nach 
Anspruch 1, gekennzeichnet dadurch, dass durch nichtsinusformige Verteilung der 
Beleuchtung des Objekts auftretende Fehler.korrigiert werden. 

8. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach 
Anspruch 7, gekennzeichnet dadurch, dass mittels Bandpassfilter hohere Haimonische 
der Grundfrequenz der periodischen Struktur aus der Helligkeitsverteilung 
herausgefiltert werden. 

9. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach. 
Anspruch 1, gekennzeichnet dadurch, dass die Verrechnung der registrierten 
Helligkeitsverteiltmgen unter Beriicksichtigung von Korrekturwerten vorzugsweise fur 
die Beriicksichtigung der Helligkeitsschwankungen der Beleuchtung, des 
Ausbleichens des Objekts" bei Fluoreszenzbeleuchtung und nichtsinusfSrmiger 
Verteilung der Beleuchtung des Objekts erfblgt, welche durch lineare Optimierung 
bestimmt werden. 

10. Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach 
Anspruch 9, gekennzeichnet dadurch, dass die lineare Optimierung auf eine Merit- 
Funktion der Fdrm 

M(&r>d;b) = * 0 \F{a} 0 \ 2 +a^ 

angewandt wird, wobei F{ a } eine Fiinktionaltransfonnierte des Vektors a ist, 9 t 

skalare Faktoren zur Charakterisierung der Beleuchtungshelligkeitsschwanktmg ist, 
d ein MaB fur das Ausbleichen des Objekts bei Fluoreszenzbeleuchtung ist, b ein 
Faktor zur Charakterisierung der nichtsinusfonnigen Verteilung der Beleuchtung des 
Objekts ist und die a { Gewichtskoeffizienten zur Anpassung an die 
Aufhahmebedingungen vorzugsweise bzgl. verschiedener Signal-Rausch Abstande 
oder bevorzugter Frequenzen sind. 

1 1 . Verfahren zur Verbesserung der Tiefendiskriminierung optischer Systeme nach 
Anspruch 1, gekennzeichnet dadurch, dass Ausbleicherscheinungen des Objekts bei 
Fluoreszenzbeleuchtung dadurch bestimmt werden, dass aus vorzugsweise 
aufeinanderfolgenden Helligkeitsregistrierungen eine lokale KoiTektiufunktion 
bestimmt wird. 
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12. Verfahrenzur Verbesserung derTiefendiskriminierung optischer System© nach 
Anspruch 1 1 , gekennzeichnet dadurch, dass diese lokale Korrekturfuriktion durch 
Mittelung uber mindestens eine Periode der periodischen Beleuchlmgsslruktur 
bestiTTimt wird. 

13. Verfahrenzur Verbesserung der Tiefendiskriminieruiig optischer Systeme nach 
Anspruch 1 1 , gekennzeichnet dadurch, dass in der Nahe von im Objekt enihaltenen 
Kanten der Axisbleichfunktion auftretendes Uberschwingungen der berechneten 

lokaleii Korrel^arfunktion durch Ersetzen mit einem Schatzwert verhindert wird. 
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